Anomalous scaling of triple correlation function of white-advected passive scalar 
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For a short-correlated Gaussian velocity field the problem of a passive scalar with the imposed 
constant gradient is considered. It is shown that the scaling of three-point correlation function is 
anomalous. In the limit of large dimension of space d anomalous exponent is calculated. 
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We consider the problem of passive scalar 9 advected 
by delta-correlated in time velocity field v [[lj . The scalar 
obeys an equation 



dtQ + vVQ = kA8 



(1) 



Here k is a molecular diffusivity, which determines be- 
havior of passive scalar for small distances. The velocity 
v is a random function to be defined below. There is no 
pumping of a passive scalar into the system, but instead 
the condition on mean value of passive scalar is imposed, 
namely it has a constant gradient in space. 



6 = gr 



(2) 



Here 9 designate the fluctuating part of a scalar field 



<0>=O. 



(3) 



By () we mean statistical average. This problem was pre- 
viously considered by Shraiman and Siggia in M , where 
they derived equation for the n-point correlation function 
and found the pair correlation function. In J3|,|| they pro- 
posed phenomenological equation describing the n-point 
correlation functions of passive scalar. 

The description of this problem, analogous to Kol- 
mogorov dimensional analysis gives correct answer for 
two-point correlation function /(ri 2 ) = (9(ri)6(r^ 
However it has a deviation from experiment data 
for higher order correlation functions. We will concen- 
trate on such deviation for the three-point correlation 
function T = (0{ri)0{r2)0{r^)) , which is the lowest that 
has anomalous behavior (i.e. different from one obtained 
by dimensional estimate) . The imposed gradient of mean 
value looks like breaking of translational and rotational 
invariance. However, the correlation functions of are 
translational invariant and the gradient of mean value 
of a scalar field will effectively play a role of anisotropic 
pumping. There is a significant difference between this 
problem and the problem with isotropic pumping. Due 
to the presence of chosen direction in space the odd order 
correlation functions do not vanish. 

We assume the velocity field to be Gaussian and delta- 
correlated in time Q. Therefore statistics of v is com- 
pletely determined by the pair-correlation function 



(v a (r 1 ,t 1 )v (r 2 ,t 2 )) = V a ^(r 12 )6(t 1 -t 2 ) 
V a(3 (r) = V 5 af3 - K a0 (r) 

Here K a @ is the eddy diffusivity tensor 



(4) 
(5) 



K af3 (r 
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2- 7 



r <* r @ 



(6) 



where d is the dimension of space and 7 is parameter 
which determines the scaling of velocity pair-correlation 
function. The parameter 7 is supposed to be between 
zero and two. Power-like behavior (j(^) is valid up to 
some r ~ L, where L is the correlation length of ve- 
locity. For r > L the correlator V al3 decreases and goes 
to zero as r — > 00. As we will see later the precise form of 
this decrease is not important. Therefore we can choose 
arbitrary form of pumping at r > L. We will believe 
K a P{r) = Vo5 a P for large r. The relation between con- 
stants is 



L 2 ^ = 



2-7 



D(d + 1 -7) 



V a 



(7) 



We assume that the cut-off length L is d-independent, so 
that Vo explicitly depends on space dimensionality. 

Our aim is to find the scaling exponent of the triple 
correlation function deep inside convective interval, that 
is for rd <C r <C L, where is diffusion scale: 



2-7 



2k(2 - 7) 
D(d-l) ' 



(8) 



The condition r r e i implies that we can disregard dif- 
fusion ||. 

From (0) one can derive the equation for the nth-order 
correlation function pip) 



£(6i...9„) =0. 



(9) 



Here the operator C contains both operators of turbulent 
and molecular diffusion: 



£ = ±£rtf(r«,)V?Vj + «X;v? 



(10) 



We will see that the scaling of the pair-correlation func- 
tion is normal that is 7. For the three-point correlation 
function the presence of zero modes of the operator C 
makes scaling anomalous that is different from 1 + 7 ob- 
tained from Kolmogorov-like estimates. The exception 
is special case when the scaling exponent of the velocity 
correlation function 7 is zero 0. Note also that for the 



1 



scaling exponent equal to two the operator L has a singu- 
larity and special treatment is needed 0. Following ||] 
we will calculate anomalous scaling exponent in the limit 
of large d. For d — > oo the problem can be solved exactly. 
Then, finding corrections in the next order over l/d one 
can calculate the anomalous exponent of the three-point 
correlator. 

To calculate the triple correlation function one should 
know the two-point correlator since it explicitly enters 
the equation. The equation for the two-point correlation 
function is 



C {p) f(r 12 ) = -C[(gr 1 )(gr 2 )} 
The operator may be written as 



(11) 



= D ^_ (r d+1 -^,.) , r d « r < L (12) 

We should match the solution obtained at r < L with 
the solution at r > L. For this region the operator has 
the following form 



(13) 



Since we will consider three-point correlation function at 
large d, we calculate pair-correlation function in the same 
limit (it can be found exactly as well). The solution of 
equation (|l|) at r < L and d 3> 1 which satisfies bound- 
ary condition at r ~ L is 



f(r) = 



2-7 
27 



r \ 7 
1/ 



1 



2 \L 



(14) 



At r > L the pair correlation function is zero in the main 
order of l/d expansion. For r deep inside the convective 
interval r 2 term is small in comparison with r 7 . There- 
fore we see that the scaling behavior of a two-point cor- 
relation function is normal, as it was shown by Shraiman 
and Siggia Note that ( p^ ) is isotropic in convective 
interval. The anisotropic part has scaling exponent larger 
then 7 and therefore in the convective interval is much 
smaller. Besides, it has additional factor l/d with respect 
to isotropic part. 

For the three-point correlation function one can derive 
the following equation 

CT = $12,3 + $13,2 + $23,1, where (15) 
$12,3 = V?/(n 2 )/ [K a ?(r 13 ) - K^(r 23 )] . (16) 

Since the preferable direction exists in the problem 
then r depends on five independent variables in any di- 
mension of space d > 2. It is convenient to choose the 
following set 



T = T (n 2 , ri3, r 23 , r 12 g, r 23 g) 



(17) 



However, this choice is is not unique because of identity 
ru + r 13 + r 23 = 0. 



Unfortunately the number of variables is too large to 
solve the equation ( |l5|) exactly at arbitrary 7 and d. 
There are several cases when the problem can be solved 
exactly: the cases of 7 = pjjojjj, of large d ||, and of 
7 = 2 [jn]]. After solving the problem exactly at these 
special values of parameters one can develop the pertur- 
bation theory for the small deviations from these values. 
In the cases of large d and 7 = 2 this perturbation theory 
can be performed. The case of 7 = corresponds to the 
strongly degenerate perturbation theory. We are going to 
solve the equation ( |l5j ) in the limit of large d, where per- 
turbation theory is regular. To do this one should keep 
in the operator (10) only terms which are the largest at 
large d and solve the equation. Then one can look for 
a correction due to a finite value of the parameter l/d. 
Justification of such a perturbation theory can be found 
in §. 

In variables (0) the main part of the operator C can 
be written as 



C 



d 2 D 
~ 7 



i<j 



"Or, 



(18) 



For large d and r -C L one can substitute ([14]) into the 
right-hand side of (Uq). The result is 



$ 



(0) 

12,3 
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{r, 2 g)r} 2 - 2 {rlr -r 2 2 r) (19) 



Thus, on the first step of iteration procedure we should 
solve the following equation 



l l o - $ 12 3 + <j> 13 2 + $23,1 



(20) 



The solution can be found by integrating over the char- 
acteristics of the operator (Eq). Important question is 
that of boundary conditions. We should supply equation 
( pp| ) with the boundary conditions at zero and infinity. 
However this equation is of the first order and only one 
condition can be satisfied. Since we omitted diffusion 
part of the operator ( |Io| ) we should pose the only condi- 
tion r — > as r — > 00. One can show || that solution 
obtained in this way can be matched with diffusion re- 
gion. 

One can write the following solution of ( po|) 



r n = 



? 2 L 2- 7 



(r 12 g) [ dt (r& + 1) 
Jo 



1-2/7 



1 d 2 

x \{t 4- r^) 2 / 7 - 1 - (t + r^) 2 / 7 - 1 + permutations (21) 

This integral diverges on upper limit and we have to regu- 
larize it. Since divergence is logarithmic then the precise 
form of regularization is not important, we just cut off 
the integral at t ~ L 1 . It corresponds to the condition 
that correlator is zero for r > L. Integral (|2l|) can not be 
calculated analytically for arbitrary 7. However we are 
interested in T in the inertial interval, that is for r <C L. 
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There main contribution to (|T]) comes from large t. Ex- 
panding integrand over r/L we find the following result 
in the main logarithmic order 



where r in the logarithm is of order r.y and 

a 1 ^ r, 



(22) 



Z 



2 7 d 



+ (n 3 £/)W 2 " r I 3 ) + ( T "239)(^ 2 - r i 3 ) 



(23) 



Note that Z is a zero mode of the operator (|lg|). It is 
necessary since logarithm enters the solution. If it were 
not zero mode, we would have logarithm at the right- 
hand side of (|16|). It is also worth mentioning that Zq 
has scaling 7 + 1 and this is the only zero mode of the 
operator C with such a scaling. 

Now we should find a correction to ( p2| ) in the next 
order over 1/d. To do this we should solve the equation 



CqTi — — £]To 



(24) 



Here C\ is the part of operator (|10|) which is propor- 
tional to d. In equation (|2j) we disregard terms which 
come from the right-hand side of the equation ( J15) s ince 
they do not contain logarithm. Solving equation ( 24j) we 
find for the main logarithmic term 



2(2 - 7 ) 2 2 (L 



(25) 



We can continue this iteration procedure and find the 
main logarithmic subsequence. As a result we obtain the 
following series 



r = ZMl+Aln(~)+^V (j 



-Z 



z„ [ - 



A 



Z , 



(26) 



where the anomalous exponent A which makes deviation 
from normal scaling is 



A = 



2(2-7) 



(27) 



We see, that the solution consists of two parts. The sec- 
ond term in the right-hand side of (^6|) has normal scaling 
7 + 1 . One can check that it is partial solution of inho- 
mogeneous equation (^). However, if it were alone, the 
solution would not satisfy boundary conditions. To en- 
sure it we have the first term in (|2g). One can see that 
this contribution gives anomalous scaling with exponent 
7 + 1 — A. This term is the solution of homogeneous 
equation CT — and therefore is the zero mode of the 
operator t §,|l]||. 

To conclude, we have shown that the three-point corre- 
lation function has the scaling exponent 7 + 1 — A which 
differs from naive dimensional estimates. This exponent 
was analytically calculated in the leading 1/d order. 
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